home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / csptri.z / csptri
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCSSSSPPPPTTTTRRRRIIII((((3333FFFF))))                                                          CCCCSSSSPPPPTTTTRRRRIIII((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CSPTRI - compute the inverse of a complex symmetric indefinite matrix A
  10.      in packed storage using the factorization A = U*D*U**T or A = L*D*L**T
  11.      computed by CSPTRF
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE CSPTRI( UPLO, N, AP, IPIV, WORK, INFO )
  15.  
  16.          CHARACTER      UPLO
  17.  
  18.          INTEGER        INFO, N
  19.  
  20.          INTEGER        IPIV( * )
  21.  
  22.          COMPLEX        AP( * ), WORK( * )
  23.  
  24. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  25.      CSPTRI computes the inverse of a complex symmetric indefinite matrix A in
  26.      packed storage using the factorization A = U*D*U**T or A = L*D*L**T
  27.      computed by CSPTRF.
  28.  
  29.  
  30. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  31.      UPLO    (input) CHARACTER*1
  32.              Specifies whether the details of the factorization are stored as
  33.              an upper or lower triangular matrix.  = 'U':  Upper triangular,
  34.              form is A = U*D*U**T;
  35.              = 'L':  Lower triangular, form is A = L*D*L**T.
  36.  
  37.      N       (input) INTEGER
  38.              The order of the matrix A.  N >= 0.
  39.  
  40.      AP      (input/output) COMPLEX array, dimension (N*(N+1)/2)
  41.              On entry, the block diagonal matrix D and the multipliers used to
  42.              obtain the factor U or L as computed by CSPTRF, stored as a
  43.              packed triangular matrix.
  44.  
  45.              On exit, if INFO = 0, the (symmetric) inverse of the original
  46.              matrix, stored as a packed triangular matrix. The j-th column of
  47.              inv(A) is stored in the array AP as follows:  if UPLO = 'U', AP(i
  48.              + (j-1)*j/2) = inv(A)(i,j) for 1<=i<=j; if UPLO = 'L', AP(i +
  49.              (j-1)*(2n-j)/2) = inv(A)(i,j) for j<=i<=n.
  50.  
  51.      IPIV    (input) INTEGER array, dimension (N)
  52.              Details of the interchanges and the block structure of D as
  53.              determined by CSPTRF.
  54.  
  55.      WORK    (workspace) COMPLEX array, dimension (N)
  56.  
  57.      INFO    (output) INTEGER
  58.              = 0: successful exit
  59.              < 0: if INFO = -i, the i-th argument had an illegal value
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCSSSSPPPPTTTTRRRRIIII((((3333FFFF))))                                                          CCCCSSSSPPPPTTTTRRRRIIII((((3333FFFF))))
  71.  
  72.  
  73.  
  74.              > 0: if INFO = i, D(i,i) = 0; the matrix is singular and its
  75.              inverse could not be computed.
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.